UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  of  THIS  page 


*TIC  FiLE  tm 


AD-A216  781 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBERS) 


6*.  NAME  OF  PERFORMING  ORGANIZATION  |6b.  OFFICE  SYMBOL 

'  If  applicable) 

PURDUE  RESEARCH  FOUNDATICN 


6c.  ADQRESS  (City,  6 late  ana  /AP  Code) 

SCHOOL  OF  AERONAUTICS  &  ASTRONAUTICS 
WEST  LAFAYETTE,  IN  47907 


PORT  DOCUMENTATION  PAGE 

ITT"  RESTRICTIVE  MARKINGS 


Ba.  NAME  OF  FUNOING/SPONSORING  8b.  OFF 

ORGANIZATION  FORCE  W«l 

OFFICE  OF  SCIENTIFIC  RESEARCH  NA 


8c.  AOORJESS  ICity.  Slate  and  ZIP  Code) 


8b.  OFFICE  SYMBOL 
If  applicable) 


3.  OlSTRIBUTIOftfl^VOlXttBiLIT.Y  .QF. REPORT  . 

distriV"  t  i 

istritutiou  unlimited. 


5.  MONITORING  ORGANIZATION  REPORT  NUM0ER(S) 

AFCHiRTR*  8  9*1  855 


7a.  NAME  OF  MONITORING  ORGANIZATION 


7b.  AOORESS  (City.  State  and  ZIP  Code) 

^Yv^>  C-  .  Q  VC' 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F49620-83-K-0008 


10.  SOURCE  OF  FUNDING  NOS. 


\6VcA  v-HO 

BOLLING  AFB  DC  10332-6448 

PROGRAM- 

PROJECT 

TASK 

WORK  UNIT 

E  LEMENT-'NO. 

NO. 

NO. 

NO. 

1 1.  TITLE  itnciuae  security  Classification*  SOLiXJ)  r^QCKLT 

SPACE  PROPULSION  STUDIES 

61102F 

2308 

'  PERSONAL  authorisi 

OBERT  L  GLICK  AND  JOHN  R  OSBORN 


TYPE  OF  REPORT 

FINAL 


16.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVERED 

from  0  1Nov82"o3  1DEC8 


44.  DATE  OF  REPORT  (Yr..  Mo..  Day) 

DECEMBER  1984 


15.  PAGE  COUNT 


17.  COSATI  CODES 


FIELO  I  GROUP  '  SUB.  GR. 


19.  SU8JECT  TERMS  ^Continue  on  reverie  if  necessary  ana  identify  by  block  number i 


19.  ABSTRACT  / Continue  on  reverse  if  necessary  and  identify  by  blocn  number i 


s 


DTIC 

ELECTEI 
JAN04 19901 

I 


b 


90  01  04  1  86 


DISTRI8UTION/A VAILABILITY  OF  ABSTRACT 

.lassified/unlimiteo  Z  same  as  rpt.  G  otic  users  □ 


21.  ABSTRACT  SECURITY  CLASSIFICATION 


22a.  NAM€  OF  R6SPONSI8LE  NOIVIOUAL 

22b.  TELEPHONE  NUMBER 

(Include  Area  Code > 

LEONARD  H  CAVENY 

(202)  767-4937 

22c.  OFFICE  SYMBOL 

AFOSR/NA 


DD  FORM  1473,  83  APR 


EOITION  OF  1  JAN  73  IS  OBSOLETE. 


UNCLASSIFIED 


SOLID  ROCKET  AND  SPACE  PROPULSION  STUDIES 


Drs.  Robert  L.  Glick 
and  John  R.  Osborn 

School  of  Aeronautics  and  Astronautics 
Purdue  University 
West  Lafayette,  IN  47907 


December  1984 


Final  Technical  Report 
Contract  F49620-83-K-0008 


Approved  for  public  release;  distribution  i_inhmitc< 


Department  of  the  Air  Force 
Air  Force  Office  of  Scientific  Research 
Bolling  Air  force  Base,  D.C.  20332 


ABSTRACT 


^"ktuaies  were  directed  at  the  deflagration  of  heterogeneous,  condensed 
media.  Relations  among  ballistic  properties  were  studied;  generalized  relations 
were  derived  and  existing  relations  for  testing  the  consistency  of  data  from 
ballistic  test  motors  were  shown  to  be  general  to  small  error;  errors  in  the 


recent  literature  were  corrected.  A  new  device  for  characterizing  the  ballistic 
properties  of  condensed  media  at  high  pressure  with  strands  was  devised  and 
explored  analytically.  3elf-pressurized  constant  pressure  operation  was  shown 
along  with  capability  to  control  the  pressure  level  with  a  simple  bang-bang  con¬ 
trol  system.  Special  configurations  to  provide  direct  little  difference  measure¬ 
ments  of  ballistic  sensitivities  were  presented.  The  Deur/Glick  serial  sandwich 
model  for  heterogeneous  propellant  combustion  was  modified  to  overcome  the 
r-7^bontinuatioiP  problem  implicit  in  that  model.  Results  showed  ignition  delays  of 
correct  magnitude  and,  with  physical  reasoning,  demonstrated  that  the  ZN 
methodology  cannot  be  applied  to  heterogeneous  propellants  in  its  present 
form.  K  o  '•/  ‘  ■  •  ,  T  ■  mV 
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1.0  INTRODUCTION 

Ccmbustion/flowfield  interactions  are  crucial  to  the  performance  of  chemi¬ 
cal  prcpulsers.  In  solid  rocket  motors  they  appear  primarily  as  condensed 
phase  media  deflagration/ ffowfield  interactions.  They  modify  mass  burning  rate, 
impact  performance,  and  dominate  combustor  stability.  Qualitative  mechanis¬ 
tic  understanding  of  the  governing  phenomena  has  not  been  achieved  to  date 
because  of  inadequate  understanding  of  the  self-induced  flowfield  and  the 
deflagration  Phenomena  of  condensed,  heterogeneous  propellant.  The  seif- 
induced  flow  field  is  unique  in  fluid  mechanics  because  deflagration  of  hetero¬ 
geneous,  condensed  media  is  spatially  and  temporally  nonsteady;  the  latter 
leads  to  the  existence  of  finite  turbulent  kinetic  energy  at  the  flow  boundary. 
Stationary  state  calculations  by  Beddini1  show  that  this  "boundary  turbulence  ' 
modifies  the  turbulence  field  and  strongly  Impacts  profile  transition  in  semi- 
enclosed  porous  ducts.  These  modifications  of  the  turbulence  field  are  of  partic¬ 
ular  importance  because  it  is  generally  conceded  that  turbulence  modification 
of  transport  phenomena  is  the  principal  mechanism  for  flowfield /deflagration 
interactions12  Thus,  the  qualitative  nature  of  the  interaction  is  clear: 
deflagration  of  heterogeneous  condensed  propellants  begets  the  flowfield  and 
supplies  boundary  turbulence;  the  flowfield  with  its  dynamics  amplifies  and 
redistributes  the  input  turbulent  kinetic  energy;  and  the  lurbuience  field 
modifies  the  deflagration  process  thereby  completing  the  loop 

It  is  clear  that  both  the  heterogeneity  of  the  condensed  phase  and 
deflagration  phenomena  at  the  scale  of  the  heterogeneities  are  important  to 
deflagration/ flowfield  interactions.  Considerable  effort  has  been  made  to 
describe  this  phenomena  theoretically.  Reasonable  success  has  been  achieved 
for  the  stationary  state  although  formulational  regions  where  the  models  fail 
miserably  exist.  Unfortunately,  little  success  has  been  achieved  with  the 
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nonstationary  problem.  Two  common  flaws  pervade  both  situations:  neglect  of 
thermal  transients  in  the  condensed  phase  at  the  scale  of  the  heterogeneities 
(observed  in  micro-thermocouple  experiments)  and  deviations  (and  nonsta- 
tionarity)  of  heterogeneity  populations  on  the  deflagrating  surface  from  the 
planar  statistics  assumed.  The  former  is  important  because  heterogeneous  pro¬ 
pellant  deflagration  is  intrinsically  nonsteady  and  at  the  scale  of  the  hetero¬ 
geneities  condensed  pnase  temperature  gradients  profoundly  effect  nonsteady 
deflagration  processes.  The  latter  is  important  because  it  means  that  what  we 
think  is  on  the  burning  surface  (planar  statistics)  probably  isn't!  Since  the  sur¬ 
face  is  rough,  one  suspects  that  the  greatest  disparities  occur  at  the  smallest 
scale  of  the  heterogeneities.  The  population  statistics  problem  also  impacts  the 
response  of  the  deflagration  rate  to  environmental  fluctuations. 

It  has  been  assumed  in  some  studies  that  the  condensed  phase  possesses 
long  range  order  so  that  deflagration  through  the  ordered  media  leads  to  an  int¬ 
rinsically  nonsteady  process  (layer  frequencies)  at  the  spatial  scale  of  that  ord¬ 
ering.  These  effects  have  been  explored  theoretically3-5  However,  since  hetero¬ 
geneous  propellants  are  prepared  by  mixing  and  significant  anisotropy  has  not 
been  observed  in  either  physical  or  ballistic  properties,  long  range  order  in  the 
condensed  phase  is  improbable.  Therefore,  the  ordered  condensed  phase  models 
are  not  particularly  realistic.  To  overcome  this  problem  Deur  and  Glick8 
created  a  serially  layered  model  with  finite  differences  in  layer  properties  and 
randomly  ordered  layers.  Calculations  with  the  model  showed  that  environmen¬ 
tal  fluctuations  did  not  correlate  with  the  randomly  ordered  structure  of  the 
condensed  phase  to  first  order  Therefore,  the  heterogeneity  response  postu¬ 
lated  by  Cohen  and  Strand5  appears  to  be  physically  implausible.  However,  the 
Deur/Glick  model  exhibited  a  fatal  flaw;  combustion  terminated  spontaneously 
with  transitions  from  fine  to  coarse  pseudo-propellant  layers  whenever  rate 
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changes  were  significant.  Therefore,  with  interesting  propeiiants  a  "continua¬ 
tion  problem"  exists. 

In  addition  to  serially  layered  models  several  versions  of  parallel  layered 
models  have  been  developed7  These  models  arise  naturally  when  one  spatially 
averages  over  the  deflagrating  surface.  Typically,  these  models  give  the  mean 
nonsteady  response  in  terms  of  the  responses  of  the  individual  pseuao- 
propellants  that  comprise  the  parallel  layers.  A  variety  of  techniques  have  been 
employed  to  deduce  the  nonsteady  behavior  of  the  individual  layers.  The  most 
general  of  the  techniques  employed  is  the  ZN  methodology3  In  this  approach 
steady  state  deflagration  information  is  employed  to  characterize  heat  feedback 
from  the  gas  phase  to  the  nonreactive  condensed  phase.  Therefore,  rate  sensi¬ 
tivities  to  pressure  and  initial  propellant  temperature  are  extremely  important 
to  these  formulations.  These  can  be  computed  from  stationary  state  models9; 
they  can  also  be  extracted  from  experimental  data  if  the  data  base  is 
sufficiently  large10 

As  noted  above,  application  of  ZN  methodology  requires  accurate  knowledge 
of  the  rate  sensitivities  n='9lnr/ dlnp]?  and  crp=tdlnr/  3T]p  .  Data  for  these  sensi¬ 
tivities  come  from  burning  rate  data  obtained  from  strands  and/or  ballistic  test 
motors.  It  is  important  to  note  that  accurate  burning  rate  data  are  necessary 
because  differentiation  (an  error  enhancing  process)  is  required.  Recent  stu¬ 
dies  by  Cohen  and  Flanigan11  have  shown  that  these  sensitivities  generally 
depend  upon  both  pressure  and  initial  propellant  temperature  ' 
n=n(p,T),<7p=.7?(p,T)  J.  Consequently,  as  previous  relations  among  ballistic  pro¬ 
perties  had  been  derived  under  a  constant  sensitivities  assumption,  new 
differential  relations  that  accommodated  variable  sensitivities  were  derived. 
Moreover,  since  burning  rate  data  are  only  available  at  discrete  pressures  and 
temperatures,  all  sensitivities  are  of  necessity  the  result  of  finite  difference 
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operations.  Therefore,  Cohen  and  Flanigan  accounted  for  finite  difference 
effects  and  found  that  significant  discrepancies  occurred  between  finite 
difference  and  differential  results.  However,  the  finite  difference  result  derived 
for  contains  a  differential  quantity  in  the  denominator.  This  suggests  that 
errors  may  exist  in  their  derivations. 

The  objectives  of  this  work  were  to  explore  the  deflagration  of  heterogene¬ 
ous  propellants  and  the  self  induced  fiowfield  in  semi-enclosed  porous  ducts  in 
order  to  increase  understanding  of  the  govermng  phenomena. 

2.0  ACCOMPLISHMENTS 

2.1  Relations  Among  Ballistic  Properties 

A  major  difficulty  with  Reference  11  is  that  confusion  has  arisen  relative  to 
differential  and  mean  values.  Cohen  and  Flanigan  take  considerable  effort  to 
illustrate  the  difference  between  these  situations  for  the  simple  case  where 
dn=dK=ciC*=dp=dTf=0  However,  in  this  simple  case  the  difference  found 
appears  to  be  due  solely  to  differences  in  the  way  ap  is  defined  fin  one  case 
<7p=Alnr/AT  v.'hile  n  the  other  crp=A r/rAT  ].  Consider  this  situation.  From 
quasi-steady  continuity 


ripAb=Atp1gc/  C* 

[l] 

and 

r^Ab=Atpzgc/  C* 

[2] 

Therefore,  for  a  constant  area  ratio  process  one  obtains 

ri/r2=pi/p2 

[3] 

Employing  the  definitions  of  rrK  and  ap  mean  values  of  these  sensitivities  can  be 
computed  by  integrating  over  the  dK-Opath  viz 


2 

7fjc=  /n-fcdT/  AT=ln(  p2/  Pi)/ AT  [4] 

i 

8 

3P= f  ffpdT/  AT 

l 


[5] 
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For  a  dK=0  patn  employ  the  identity  jp=tfK-n^K  and  obtain 


<7p=[in(r2/  ri)-nln(p2/  Pi)]/AT 


With  [3,4.6] 


Since  this  is  precisely  the  result  obtained  for  a  differential  change  under  these 
assumptions,  difference  and  differential  expressions  are  identical  for  this  res¬ 
tricted  situation  in  direct  contradiction  to  Reference  11! 

For  the  relation  among  ballistic  properties  in  ballistic  test  motors  it  can  be 
shown  in  general  that 

ffp+ffc+'dln(pB-p,WTjK 

„K= - -  ,3] 

To  obtain  the  Cohen/Flanigan  expression  expand  n  in  a  Tavior  series  about  the 
reference  state  0°  ,  truncate  that  expression  at  first  order,  and  employ  the 
identity 


[an/  3T]K=Lan/  9T]p-7rK[6n/  ainpj-r 


to  obtain 


gp-j-n-c+[ain(pc-p;)/aT;K _ 

K  l-n0-[(an/aT)p+rrK(an/ainp)-]0[T-T°]  ’■*  J 

If  the  finite  temperature  difference  is  replaced  by  a  differential  temperature 

difference,  the  Cohen/ Flanigan  result  is  recovered.  However,  this  is  still  a 

differential  expression.  That  is.  the  rrK  computed  with  [10]  is  not  a  mean  value 
for  the  T°  to  T  range. 

To  obtain  an  appropriate  mean  value  for  ttk  the  differential  expression  from 
[8]  must  be  integrated  over  the  path  of  the  process.  Consequently, 

1-n  L  1 

If  n,  C*.  etc  are  known  functions  of  p,T  this  differential  expression  can  be 

integrated  over  the  range  from  state  1  to  state  Z  and  the  mean  value  of  tt* 
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determmed  from  the  definition.  This  methodology  allows  for  variable  sensitivi¬ 
ties  in  a  natural  way. 

In  a  recent  work  Gaunce  and  Osborn12  derived  another  expression  for 
rrK=func(ap  etc)  that  included  sensitivities  of  sensitivities.  This  expression  was 
reported  to  be  the  only  completely  general  expression  in  existence.  However,  it 
was  found  to  be  an  algebraic  identity  of  the  expression  derived  previously  by 
Reference  13.  From  St.  Robert's  law  lnr=lnc+ninp  so  that  differentiation  gives 

<7K=[3lnc/  9T]K+n[5lnp/  3TjK+lnp.3n/ 3T]k  [12] 

With  the  general  relation  aK=ap+nrrK  [12]  becomes 

[31nc/  3T][<=crp-lnp[3n/  3T]k  [  13] 

Substituting  this  result  into  the  expression  of  Gaunce  and  Osborn  yields  [8], 

Therefore,  the  results  are  mathematically  identical!  Unfortunately,  they  are  not 
identical  at  a  practical  level  because  sensitivities  of  sensitivities  introduce  addi¬ 
tional  error  when  operating  with  real  ballistic  data  that  is  of  necessity  imper¬ 
fect.  This  is  easily  demonstrated  by  computing  rrK  from  data  supplied  by  Gaunce 
and  Osborn  note  that  the  absence  of  error  in  Table  5  for  the  Gaunce  and  Osborn 
expression  results  from  the  arbitrary  definition  of  null  error  for  that  result!] 
which  gives  different  numerical  results  for  mathematical  identities.  Conse¬ 
quently,  Gaunce  and  Osborn’s  expression  for  ttk  is  overly  complex. 

The  necessity  for  accurate  and  consistent  ballistic  data  for  sensitivity  cal¬ 
culations  suggests  that  the  procedure  developed  and  employed  by  Geckler  and 
Sprenger  be  used  to  improve  the  quality  of  the  data  base  by  eliminating  incon¬ 
sistent  data.  It  has  been  shown  in  Reference  13  that  this  methodology,  derived 
for  constant  sensitivities,  is  general  and  applies  to  variable  sensitivities.  For 
details  see  Reference  13. 


2.2  Ballistic  Test.  Devices 
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As  noted  previously  burning  rates  are  determined  experimentally  in  either 
strand  burners  or  ballistic  test  motors.  Strand  burners  connected  to  large 
surge  tanks  yield  data  at  essentially  constant  pressure.  However,  operation  of 
this  type  device  at  high  pressures  is  difficult  if  not  impossible  because  of  purge 
difficulties.  Therefore,  strand  data  are  collected  at  high  pressure  in  closed 
burners  that  of  necessity  do  not  maintain  constant  pressure,  '"'his  is  not  a  prob¬ 
lem  with  ballistic  test  motors  because  they  are  self  pressurizing  as  long  as  noz¬ 
zle  plugging  is  not  a  problem.  However,  in  all  cases  only  absolute  rate  data  are 
obtained  so  that  differences  required  to  obtain  the  sensitivities  introduce  appre¬ 
ciable  error.  A  new  test  device  -  a  hydraulic  strand  burner  -  has  been  devised 
that  permits  self  pressurized  strand  burning  at  elevated  pressures.  Moreover,  a 
manifolding  arrangement  has  been  devised  that  permits  direct  measurement  of 
the  little  difference  in  burning  rate  at  different  states  appropriate  for  either  r 
or  <7p  measurement.  The  concepts  are  well  described  in  Reference  14 

2.3  Combustion  Modeling 

The  serially  layered  model  of  Deur  and  Glick  was  modified  to  overcome  the 
"continuation"  problem  noted  above.  The  modifications  were  limited  to  bimodal 
propellants  although  the  methodology  could  be  applied  to  the  polvmodal  situa¬ 
tion.  In  a  bimodal  formulation  of  coarse  and  fine  particulates  the  size  dispartitv 
leads  to  a  significant  probability  that  both  coarse  and  fine  particulates  will  be 
adjacent  to  fines.  As  ignition  delay  increases  with  particle  size,  the  fine  pseudo¬ 
propellant  will,  to  a  reasonable  approximation,  behave  like  a  homogeneous  pro¬ 
pellant.  Therefore,  the  coarse  pseudo-propellant  is  literally  surrounded  by  an 
ocean  of  fine  pseudo-propellant  that  behaves  like  homogeneous  propellant.  Pre¬ 
vious  calculations0  indicated  that  spontaneous  extinguishment  occurs  only  in 
transitions  from  fine  to  coarse  pseudo-propellant.  Therefore,  when  the  coarse 
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pseudc-propeiiani  extinguishes,  the  fine  pseudo-propellant  continues  burning. 
Consequently,  the  extinguished  coarse  particle  sees  a  surrounding  flow  from  the 
fine  pseudo-propellant  and  heat  transfer  from  that  flow  will  add  energy  to  the 
unignited  particle.  The  flow  situation  will  be  similar  to  separated  flow  at  the 
trailing  edge  of  a  sphere. 

This  physical  situation  was  modeled  and  embedded  in  the  basic  Deur/Glick 
model.  Calculations  showed  that  oxidizer  particle  igmtion  delays  computed  in 
this  manner  were  of  proper  magnitude15  Unfortunately,  there  was  insufficient 
computer  time  for  the  random  process  to  achieve  stable  means  for  the  pressure 
coupled  response  function.  However,  a  reasonable  estimate  was  obtained  by 
smoothing  the  computed  results.  This  showed  that  multiple  maxima  type 
response  was  not  achieved.  This  is  not  a  complete  surprise  if  one  considers  the 
physics  of  the  process.  In  this  simple  model  the  fine  pseudo-propellant  behaves 
like  a  homogeneous  propellant  -  it  is  always  burning !  However,  the  coarse 
pseudo-propellant  exists  in  two  states;  unignited  and  ignited.  In  the  unignited 
state  it  is  dormant  and  does  not  contribute  to  the  pressure  coupled  response. 
However,  the  ignited  state  is  not  an  equilibrium  state  because  of  the  augmented 
energy  store  m  the  condensed  phase  accumulated  during  the  igmtion  delay. 
This  excess  of  sub-surface  energy  will  force  the  igmted  state  burning  rate  to 
higher  than  equiubrium  values.  Consequently ,  the  coarse  pseudo-propellant  will 
behave  as  an  inert  diluent  during  its  ignition  delay  period  and  somewhat  Like 
fine  pseudo-prope llant  during  its  ignited  state.  Note  that  the  balance  between 
these  processes  will  depend  upon  particle  size,  pressure,  initial  temperature, 
and  the  deflagration  characteristics  of  the  fine  pseudo-propellant. 

The  realization  that  heterogeneous  propellant  deflagration  is  nonsteady  at 
the  single  particle  level  has  considerable  importance.  For  example,  it  implies 
that  the  Z N  memodology  in  its  present  incarnation  is  not  applicable  to  hetero- 
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geneous  propeiiants.  To  see  this  note  that  the  sensitivities  employed  with  the  ZN' 
method  are  those  of  the  mean  state.  That  is,  the  sensitivities  one  obtains  from 
experimented  data  are  of  necessity  those  for  a  mean  of  ignited  and  unignited 
states.  However,  only  the  ignited  states  have  reality  for  nonsteadv  response! 
Therefore,  the  present  ZN  methodology  is  not  applicable  to  parallel  layer 
models. 

2.4  Miscellaneous 

Rocket  motors  were  designed  for  concurrent  mean  and  nonsteady  pressure 
and  nonsteadv  velocity  measurement.  These  designs  were  implimented.  Unfor¬ 
tunately,  reduced  data  were  not  achieved  because  of  a  data  acquisition  system 
failure;  power  conditioning  at  the  T3PC  is  marginal  at  best;  power  "glitches'  are 
not  conducive  to  microcomputer  operation. 

3.0  CONCLUSIONS 

The  worxs  performed  in  this  study  have  shown  that  the  finite  difference 
/TK=-Kiap,etc)  relation  derived  by  Reference  11  is  in  actuality  a  differential  rela¬ 
tion;  a  methodology  for  obtaining  the  correct  relationship  was  derived.  Rela¬ 
tions  among  ballistic  properties  were  derived  and  it  was  shown  that  the  con¬ 
sistency  procedure  devised  by  Geckler  and  Sprenger  for  constant  sensitivities 
applies  to  the  variable  sensitivity  situation.  A  new  experimental  device  for  self- 
pressurized  deflagration  at  either  constant  or  controllable  pressure  was  devised 
and  explored  analytically.  This  device  has  capabilities  for  measuring  perfor¬ 
mance  related  parameters  and,  in  special  configurations,  is  capable  of  direct 
measurements  of  burning  rate  differences.  The  serially  layered  model  of  Deur 
and  Ghck  was  modified  to  overcome  its  ’continuation"  problem.  Results  and 
physical  reasoning  indicate  that  the  present  ZN  methodology  is  not  applicable  to 
heterogeneous  propellants. 
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